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Abstract- Aim of the paper is the qualitative analysis of a quasi-linear parabolic 
third order equation, which describes the evolution in a large class of dissipative 
^ • models. As examples of some typical boundary problems, both Dirichlet's and 

Neumann's type boundary conditions are examined. In the linear case, the related 
Green functions are explicitly determined, together with rigorous estimates of their 
tN ■ behavior when the parameter of dissipation e is vanishing. These results are basic 

Cn . to study the integral equations to which the non linear problems can be reduced. 

^^ I Moreover, boundary layer estimates can be determined too. 

Riassunto- Oggetto del lavoro e' I'analisi qualitativa di un'equazione parabolica 
p\^ , quasi lineare del terzo ordine, che descrive I'interazione tra propagazione ondosa 

^ I e diffusione in ample classi di modelli dissipativi. Quali esempi di problemi al 

contorno, vengono prese in esame condizioni sia di tipo Dirichlet che di tipo Neu- 
mann. Nel caso lineare, si determinano esplicitamente le relative funzioni di Green, 
insieme ad alcune stime rigorose del loro comportamento quando il parametro e 
della dissipazione tende a zero. Queste stime costituiscono la premessa per I'analisi 
delle equazioni integrali cui e' possibile ridurre i problemi al contorno non lineari. 
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Introduction 



A great deal of models of applied sciences are described by the parabolic 
equation: 



(1.1) CsU = euxxt + c^Uxx - utt - 2aut = -/. 

The constants a, €^,8 are all positive and they assume various meanings 
according to physical problems. As for the source /, it can be linear or not. 

For instance, the equation (1.1) is involved in the generalized Maxwell- 
Cattaneo system of equations [l]-[3], in problems of viscoelastic media of 
Kelvin- Voigt type [4] , or for the study of solids at very low temperatures [5] . 
Further applications arise in the study of viscoelastic plates with memory, 
when the relaxation function is given by exponential functions. ([6] and 
references therein). 

A typical example of the non linear case is the perturbed sine- Gordon 
equation which models the flux dynamics in Josephson junctions in super- 
conductivity [7] [8]. In this case, the terms euxxt and aut characterize the 
dissipative normal electron current flow along and across the junction. 

As for the practical applications of superconductors, many areas are 
involved. In medicine, for instance. Magnetic Resonance Imaging (MRI) 
has been used since 1977 and is still improving [9]. Referring to the electric 
power systems, the high temperature superconductor cables are likely to 
lead to a lot of benefits as regards the current carrying capacity and for 
reducing electrical losses. [10] [11]. 

As for typical boundary value problems related to the equation (1.1), 
both the Dirichlet conditions and Neumann conditions have interest for 
practical applications. For instance, in superconductivity, the first case 
can be referred to periodic conditions according to annular geometry of 
junction [12] [13], while in the other case, the phase gradient, proportional 
to the magnetic field, is specified. ([14]-[17]). When the source term / is 
linear, all these problems can be explicitly solved by means of the Fourier 
method. The solutions are determined in sect. 2-3, together with the related 
Green functions Ge,-fCe. 

When the function / is non linear, then Gs , i^e represent the kernels 
of the integral equations to which the above mentioned boundary value 



problems can be reduced. For this, a rigorous analysis of the behavior of 
these kernels when e — )■ and t — t- oo is achieved in sect 4. At last, as 
first application, the influence of the dissipation on the wave behavior is 
estimated by an asymptotic approximation uniformly valid also for large t 
(sect. 5). 



2 Statement of the problem 



If Us{x, t) is a function defined in the strip 

n = {{x,t) : < z < TT, t > 0}, 

let Ve the initial- boundary value problem related to equation (1.1) with 
conditions 

(2.1) Ue{x,0) = fo{x), dtUe{x,0) = /i(x), X G [0,7r], 

(2.2) Ue{0,t)=^{t), Ue{TT,t)=ll,{t), i > 0, 

where fo,fi,ip,^ are arbitrary date. 

The boundary conditions (2.2) represent only an example of the anal- 
ysis we are going to apply. Equally, flux-boundary conditions or mixed- 
boundary conditions can be considered too. So, another example is given 
by the problem Ti^ defined in Q by (1.1)-(2.1) together with the Neumann 
conditions 

(2.3) a,n,(0,t)=(^i(t), a,n,(7r,t) = Vi(t), t>0. 

When e = 0, the parabolic equation (1.1) turns into the hyperbolic 
telegraph equation 

(2.4) CqUo = {c^dxx -dtt -2adt)uo = -f{x,t,uo) 



and the problem Ve changes into a problem Vq for uo{x,t) which has the 
same initial-boundary conditions (2.1) - (2.2) of Vs. When the source term 
/ of (2.4) is linear (/ = f{x,t)), Vo is explicitly solved by means of the 
well-known Green function: 



(2.5)Go(x,e,t) 



2 _„ f ^ sin (i \/ (?rfi — a^ 



— e 



oo 
n=\ 



yjc^v? — a? 



sin(nx) sin(n^) 



In order to estimate the influence of the dissipative term e Uxxt on the 
wave behavior of mq ■, the difference 

(2.6) v{x,t) = Ue -uo, 

is to be evaluated and so the following problem A must be analyzed 

e Vxxt + c^ Vxx - vtt - 2avt = -F{x, t, uq, v) (x, t) £ il. 

(2.7) < v{x,0)=0, vt{x,0)=0, xG[0,7r], 



[ viO,t) = 0, vi7r,t) = 0, 



t> 0. 



The source term F is given by 

(2.8) F = f{x, t, Uo + v) - f{x, t,uo) + e uqxxu 

while, in the linear case, it is f = f and F = euQxxt- 

As for the problem ^g, instead of (2.7)3 the following conditions 



(2.9) t;^(0,t) = 0, f^(7r,t) = 0, t > 0, 

must be specified. 



3 Linear case and explicit solutions 



Let z{s) the Laplace-trasform of the function z{t) and let 



(3.1) a{s) = 



es + c^ 



When F is Unear and the Laplace trasform is applied to the problem A, 
the transform v{x,s) of the solution v{x,t) is given by 

(3.2) v{x,s)= rGe{x,ts)F{ts)dt 

Jo 



where 



(3.3) 4= ^.J^^2^ [9{\x-^\,a)-g{\x + ^\,a)] 



and 

/ ,N ^/ N cosh [(vr — y) a] 

3.4 9{y,c7) = ^^ ,. % \ 

a senn[TTa) 

But, for y G [0,27r], it results [18]: 



n=l 

and COS (n|x — .^1) — cos(n|x + ^|) = 2 sin (nx) sin (n^). So, by (3.3)- (3.5), 
it follows 

(3.6) G.(x,e,.) = - E «-«)«-(--) 



vr -^ — ' s^ + 2as + (es + c^)?^^ ' 

n=l ^ ' 



which represents the L - transform of the Green function related to problem 
A. 

By means of elementary formulae, one deduces 



2 °° 



^ 1 

n=l 



(3.7) Geix,^,t) = - ^ hnit,e) sin(nx) sin(n^), 

with 

(3.8) hnit,e) =e-(§"'+'^)* 



£„24-.u sin[tw'c2n2-(|n2 + a)2; 



Jc^n^ — (|n2 + a)2 
Then, the esplicit solution v of the problem A is: 

(3.9) v{x,t)= f dT r F{^,t) Ge{x,^,t-T)d^ 

Jo Jo 

with the Green function G^ defined by (3. 7), (3. 8). 

The formal analysis developed so far can be justified as follows. Refer- 
ring to (3.9)- (3.7), the terms 

(3.10) Fn{t) = - r F{tt) sin {nOd^ 

vr Jo 

represent the Fourier coefficients of the sine series of the function F{x, t): 

oo 

(3.11) F{x,t) = J2 Pn{t) sin(nx) 

n=l 

and the rapidity of pointwise convergence of this series depends, of course, 
on the properties of the source F. For instance, it can be sufficiently as- 
sumed that F, Fx, Fxx are continuous in (0, vr) and more 

(3.12) F(0,t) =F(7r,t) = 0. 



Then, the convergence of (3.11) is uniform everywhere in [0,7r] and, 
further, it results: 

(3.13) F„(t) = -^2 - r%(e,t) sin«) d^. 

n"= vr Jo 

As consequence, if one puts: 

(3.14) Vnit) = hn* Fn= f F„(r) hnit-T) dr, 

Jo 

the solution (3. 9), (3. 7) represents the Fourier sine expansion of v{x,t): 

(3.15) v{x,t) = 22 '^n{t) sin(nx). 

n=l 

Theorem 3.1- When F{x,-) G C^ (A) and satisfies 3.12, the solution 
v{x,t) of the problem A can be given the form : 

(3.16) v{x, t) = - f dr r % (e, r) H, (x, C, t - r) d^, 

Jo Jo 

where Hp is 



(3.17) H,{C,x,t) = - y IhL^lll sin(nx) sin(ne) 

vr -^ n^ 

n=l 

and the convergence of the series is uniform everywhere in [0,it]. m 

Remark 3.1- Theorem 3.1 can be applied also to the problem ^g, pro- 
vided that the Green function G^ is substituted by the following function: 



h (t) 2 °° 

(3.18) i^e(^,x,t)= + — / hn(t,e) cos(nx) cos(n^), 

vr vr "^ 

n=l 



4 Estimates and properties of the series H^ 



The arguments of sine functions in (3.8) are real parameters when Ki < 
n < K2, with 



c,, L 2ae , ^, c ,^ I 2ae 



(4.1) i^i= -(l-Wl-_- ), K2=-(l + \l 



e V c^ eve"' 



So, if a < c and N = [K2] , the /i„ 's in (3.8) contain trigonometric 
functions for 1 < n < iV and hyperbohc functions for n > A^ + 1. Oth- 
erwise, if a > c, the trigonometric case is related only to [Ki] < n < N. 
This distinction is unimportant to what we are going to demonstrate; how- 
ever it holds also for the Green function Go defined by (2.5) and related to 
the problem Vq. 

Let gn{x,^) = {2/tt) sin(nx) sin(n^) and consider the series 



°° sin ( t \J (P'v?' — a^ 



I , i^\ r-r —ni \~^ ''^^i I (, V <-- 't LI J , ^, 

(4.2) Ho = e^' Y. A. .2^2 J 9nix,0 



n=l 



v~ 



n^ yc^n^ — a 



deduced from iJg setting formally e = 0. 

In order to estimate the difference H^ — Hq, let 



^ 2. d2 _ „2^2 a2. 1,2 _ „2„2 „2 



(4.3) An =a + -n^; Bf^ = c^n^ - Af^; h^ = c^n^ - a 

and 



(4.4) r„ 



hn{t,e) hn{t,0) _ _At sin{Bnt) _^^ sin{bnt) 



'n? n^ n^ Bn n'^ bn 

with hn{t,e) defined by (3.8). It results: 



N oo 

(4.5) He-Ho =Rl +R2 =Y1 rn{t,e) gn + Yl ^nit,e) Qn. 

n=l n=N+l 

If cq denotes the Euler constant (cq — 0, 5773) and k an arbitrary 
constant such that < /c < 1 , let h\ = c^ — a^ and 

4.6 p{t) = -r-{2 + at), c,{e) = .: '. • 

Ol TT C(l — K) 

Further, let C2 = (1/3) r7T,ax(l, vr/cfei). Then one has: 

Lemma 4.1 - For all t > 0, a; G [0,'7r], when e is vanishing, the 
following estimates hold: 



(4.7) I ^1 I < £ ci(e) /9(t)e' 



-at 



(4.8) \R2\< ec2 [e"'^* + 06"-^], 

where 6 denotes the fast time tje and e ci (e) vanishes with arbitrary order 
k <1. 

Proof- Referring to the trigonometric terms related to i?i, defined in 
(4.5), by means of the Laplace transform, by (4.4) one deduces that 



(4.9) r„(s,e)=--^ 



bn (s + a)2 + bl (s + ^„)2 + Bl 
hence 



(4.10) r„(i,e) = -l[e-«*sin(&„t)] * [e-^"*(^ sin(5„ t) - cos(i?„t))]. 
When this convolution is made explicit, by elementary estimates one has 



^ -at 



(4.11) \rn{t,e)\ < -/9e-"* n G [1, iV] 

n 

and so 



2 ^ f 2p- 1 

(4.12) p-^ e-' \ Ri\ < - J2 - < — {CO + — + InN). 

IT , n TT ZiV 
n=l 

For each positive constant /? , one has InN < /3 N^'^ so that for /3 = 

k^^ {k < 1) the estimate (4.7) follows, with p and ci defined by (4.6). As 
for i?2 one has: 



(4.13) i?2= E ^/^n(t,e)<7n-e-'^* £ ^^^ff„ = i?^ - i?^' 

n=N+l'^ n=N+l ^ "" 

where the terms /i„ defined in (3.8) represent now hyperbolic functions 
{B'l < 0). For this it results: 



(4.14) Kit,e) < t e J^I^TS^ Vn > A^ + 1 

and A^ + 1 > 2c /e for e < 2(c — a). As consequence: 



2 t _ c2 t ^ 1 



(4.15) I -R2 I < — e" - - y -^ = (1/3) e e~ - 

n=l 

with 6* = t/e. At last, as 6„ > 61 n, one deduces that 



(4.16) e-|i?^'|< :;!: E i ^^(^ + 1)-^^ "' 

n=Af+l 



^&i ^rr;li « 36i 3c 61 



and (4.15)-(4.16) imply the estimate (4.8). 
Referring to (4.8), let observe that 

10 



(4.17) ee-^^ < {4/c^e) e"^^ 

and let b = ?7im (a, c^/4e).Thus, by (4. 6), (4. 7), (4. 8) one has 

(4.18) \ Ri \ + \ R2 \< e''r{t)e-^* 
with 

(4.19) r{t) =e^'^ [ci{e) p{t) + 02(1 + 4/0^6)]. 

Then, the following theorem can be stated. 

Theorem 4.1 - Whatever the positive constant k < 1 may he, for 

all t > and x G [0,7r], it results: 

(4.20) \ He - Hq\< ^e^{l + t + t^)e~^\ 

where the constant 7 depends only on k, a, c. 

Remark 4.1- The asymptotic analysis of this section and the results 
of theorem 4.1 can be applied also to the function i^^ defined by (3.18) and 
related to the problem Tie- 

5 Conclusions 



To outline a first application and to avoid too many formulae, let con- 
sider only the term depending on the source. Then, referring to the prob- 
lems Ve and Vo and putting AH = H^ — Hq , it results 

(5.1) ue- uo =- r di f /g5(e,r) AH{x,C,t-T)dT. 

Jo Jo 



11 



By assuming that / G C^(il) and that fxx is bounded also when 
t — )■ oo , let 

(5.2) ||tt/|| = sup \fxxix,t) \. 

n 

As consequence of theorem 4.1, when e — )■ , the following rigorous ap- 
proximation holds: 



(5.3) Ue = uo + e'^ r y{x,t) G n 
where the error r is such that 

(5.4) \r\ < 71 \\uf\\ Vi> 

and the constant 71 depends only by a, c, k. So, the error of the approxi- 
mation is neglegible also for large t (t — )■ 00). 

When / is non linear, an integral equation like 



(5.5) V = r d^j AH{x,C,t - T)F[(,T,u{(,T),Ur{C,r),Ui:{C,r)]dT 
Jo Jo 

must be analyzed. These applications will be dealt successively. 
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